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SUMMARY 

It is shown in this paper that a conventional electronic frequency com- 
pensator does not provide adequate compensation near the resonant frequency 
of a lightly damped second order system, such as the drag force anemometer. 
The reason for this is discussed, and a simple circuit modification is pre- 
sented which overcomes the difficulty. The improvement is shown in theore- 
tical frequency response curves as well as in the experimental results from 
some typical drag force anemometers. 

INTRODUCTION 

The use of electrical or electronic means to extend the frequency re- 
sponse of a transducer is not new. For instance, passive electrical net- 
works have long been used to improve the response of thermocouples (ref. 1), 
and a frequency compensation is a practical necessity in hot wire anemometry, 
raising the useful frequency range from a few hundred Hz to nearly 100 kHz. 

A similar situation has arisen in the application of the drag force ane- 
mometer (ref. 2), which is a probe that is used to measure either velocity 
head or flow angle in moving fluids. Because of the design of the device, 
which is described in greater detail in the next section, it has inherently 
2ood frequency response. Some designs have a natural frequency as high as 
43 kHz, but there arise situations such as flow monitoring in jet engines, 
where frequencies as high as 100 kHz are of interest. This paper describes 
the problems encountered in applying a frequency compensator to the drag 
force anemometer, describes the remedy, and describes the results of using a 
modified compensator on sample drag force anemometers. 

By way of background, the paper begins with a brief description of the 
drag force anemometer and then a brief discussion of second order systems in 
general. Next, electronic frequency compensation is discussed along with 
the problems encountered in applying compensation to the drag force anemome- 
ter. It is then shown that the problem may be overcome by means of a simple 
change in the compensator circuit. Both calculated and experimental results 
are presented which show the improvement. The detailed mathematics are 
given in an appendix. 

DESCRIPTION OF THE DRAG FORCE ANEMOMETER 

The drag force anemometer has found application in many different flow 
measurement areas. Some of the applications are described in reference 2. 







The drag force anemometer is a cantilever beam, with strain gauges 
mounted near the base. When the free end of the beam is exposed to a moving 
fluid, the fluid exerts a drag force F(t) on the beam. If the beam orien- 
tation is approximately normal to the flow, F(t) is proportional to the 
velocity head of the moving fluid. If the beam is oriented so that the flow 
is roughly parallel to the beam, F(t) is proportional to the product of the 
velocity head and flow angle (see figure 1). Thus, in one probe, one has a 
tool for measuring both velocity head and flow angle. 

The drag force exerted by the fluid on the beam causes a strain in the 
beam and a resistance change in the strain gauges. This is turned into an 
electrical signal x(t) by the associated electronics. 

It is shown in reference 2 that the dynamic response of the anemometer 
is described by the equation. 


x/o)^ + 2 c + X = K F(t) (1) 

The parameter f^ = is the natural frequency, c is the damping 

coefficient, and k is a proportionality constant. 

If fp is large compared to any frequency in the flow, the first two 
terms in equation (1) may be omitted, and equation (1) reduces to x = <F(t); 
that is, x(t) is directly proportional to F(t) and the output is an accu- 
rate representation of the dynamic velocity head (or product of velocity head 
and flow angle, as explained above). Hence, one way to improve the frequency 
response is to increase the natural frequency. The natural frequency can be 
increased by making the beam of a high specific stiffness material such as 
silicon or alumina, instead of a material like stainless steel, whose speci- 
fic stiffness is lower. The natural frequency can also be increased by mak- 
ing the beam shorter and thicker, but this reduces sensitivity. 

The drag force anemometer shown in figure 2 has a beam made of silicon 
and a natural frequency of 42.8 kHz. This is the highest natural frequency 
of any of the beams built at the Lewis Research Center. 

'' DYNAMIC BEHAVIOR OF SECOND ORDER SYSTEMS 

The frequency response of the drag force anemometer and other second 
order systems (i.e., systems described by linear second order differential 
equations such as eq. (1)) may be examined by letting the driving function 
icF(t) =^o with J’q constant. With this function on the right side 

of equation (1), the solution x(t) is 


x(t) = A sin lot + B cos ut = R(io) sin [wt - <p(to)] 


R(u) , tan [ 9 (u)] = B/A 
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The amplitude R and phase shift <p are 



The amplitude response defined as 20 log R/Sq and q> are plotted in 
figures 3 and 4 for the case f^ = 1300 Hz and c = 0.024. These values are 
for a drag force anemometer similar to the one in figure 2, except that the 
beam is slightly longer and made of stainless steel rather than silicon. 

The figures for different damping coefficients would be essentially the 
same, with the peak near f^ being higher and sharper for smaller values 
of c; such figures are widely available throughout the technical literature. 

These figures show that the amplitude response and 9 are strongly 
dependent on frequency above f/f^ = 0.2. In fact, for this example, the am- 
plitude increase is about 4 percent (0.4 dB) at f = 0.2 f^ and is about 10 
percent (0.8 dB) at f = 0.3 f^. If no correction is made, these changes 
become measurement errors, in addition to any other system errors present. 
They are an inherent characteristic of second order systems. 

Even if higher frequencies are not of interest in the measurement, there 
may be frequency components in the flow that can excite a drag force anemo- 
meter beam to resonance. The result is a signal in which the low frequency 
information is obscured by large amplitude oscillations at the natural 
frequency. 

COMPENSATION: PRINCIPLES OF OPERATION 

Since the problem with using the drag force anemometer to make dynamic 
flow measurements is that the amplitude response and phase vary with fre- 
quency, a solution would be to feed the signal into a device (amplifier) 
whose gain and phase shift varied with frequency in a manner exactly the 
inverse of the way they vary in the anemometer. The large increase in ampli- 
tude response near the natural frequency would be attenuated, but for higher 
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frequencies the gain would increase, making up for the decreasing sensitivity 
of the anemometer. 

This is what happens in the compensator, but instead of viewing it as a 
device whose gain varies with frequency, it is easier to understand the prin- 
ciple of operation by considering what it does to the signal it gets from 
the anemometer. As is shown in figure 5, the anemometpr output x(t) is 
fed into the compensator which electronically performs’ on x(t) all of 
the mathematical operations indicated^in equation (1). The ideal differen- 
tiators introduce a phase shift^of 90* at the output of the first differen- 
tiator and a phase shift of 180* at the output of the second differentiator, 
both referred to the input. This compensation completely eliminates the 
peak centered at f = fn shown in figure 3. The output is icF(t), which 
is the desired signal. This topic is discussed more fully in reference 3. 

COMPENSATION: SOME PRACTICAL CONSIDERATIONS 

Since real circuits and real signals are not ideal but always have noise 
associated with them, some practical considerations arise in the application 
of compensators. One consideration is that it is not possible to extend the 
frequency range indefinitely. As the sensitivity of the anemometer continues 
to drop, a point will be reached where the signal from the anemometer will be 
buried in the noise at that frequency and the compensator will not distin- 
guish between them. Another consideration is that because of the differen- 
tiators it contains, a compensator tends to act as a noise source in its own 
right. As is explained in reference 4, the gain of a differentiator in- 
creases linearly with frequency, so it amplifies any high-frequency noise. 

A standard operational amplifier differentiator circuit is shown in fig- 
ure 6(a). In order to reduce the noise, the circuit must be modified so 
that the gain does not continue to increase with frequency. The usual modi- 
fication is to add a small resistor in series with the input capacitor, and 
to add a small feedback capacitor. The modified circuit is shown in figure 
6(b). These changes reduce the gain at high frequencies, but permit approx- 
imate differentiation at lower frequencies. All the compensators discussed 
in this report, whether called "conventional" or "modified", are assumed to 
be built with modified differentiators. 

A compensator built using the modified differentiators shown in figure 
6(b) had a peak-to-peak noise of about 5 mV (0 to 100 kHz) with the input 
shorted. This is acceptable in most applications of the drag force anemo- 
meter. 


ADDITIONAL CONSIDERATIONS FOR SYSTEMS WITH LOW DAMPING 

Because of the modifications required in the differentiators, their out- 
put is no longer exactly proportional to the derivative of the input; there 
is now a small error. When these modified differentiators are incorporated 
into a compensator, there is also an overall compensator error. The error 
is most important for frequencies near f^,, and the size of the error is 
related to the size of the damping coefficient. 

The example plotted in figure 7 shows that if the damping coefficient is 
not too small, the conventional compensator of figure 5 is able to achieve 
adequate compensation by maintaining the amplitude response flat to within 
1 dB for frequencies out to about two and one-half times fp. Note that 
unlike the results obtained with perfect differentiators, the peak is not 
completely eliminated since cancellation occurs only at fp and not over 
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the entire frequency range. Figure 8, however, shows that for a system with 
a smaller damping coefficient the conventional circuit (with modified dif- 
ferentiators) is unable to provide adequate compensation due primarily to an 
inability to achieve cancellation at f = fp. 

Fortunately, it turns out that just by removing the inverter from the 
circuit of figure 5, the cancellation at fp can be restored. The details 
are in appendix B, but the idea can be discussed briefly here. 

To see what causes the problem, suppose that the compensator input is 
ei = j/sin u)t. Let the output of the first differentiator be e2 and the 
output of the second be e3, and let the damping adjustment be set at 2k. 

With a compensator input ej = ji/ sin wt, equation (610) of appendix B 
gives for ez 


e, = r s/— cos (cot - o) = — sin 


and equation (Bll) of appendix B gives for 03 


cot 




(4) 


2 2 

e^ “ sin (cot - 2o) = sin [cot - (it + 2e)j 

n n 


(5) 


A perfect differentiator introduces a phase shift of n/2 between its input and 
output, so that 0 is a small additional phase shift due to the noise reduction 
circuitry. For a perfect differentiator, e = 0. 

The output of the compensator of figure 5 is ej = e^ - 2 ke2 +63. 

Since o is small, cos o = 1, sin o = e, and 


) sin cot + 2k — cos <ot 
0) 

2 

+ 2 k o — j/ sin cot + 2 e cos cot (6) 


Without the noise reduction modifications on the differentiators, the last 
two terms in equation (6) are not present, since e = 0. At f = fp, the 
first term vanishes, and if k = c, the damping coefficient of the anemo- 
meter, ej; = 2c ja/’ cos cop t. This is what the compensator output should be 
for an input aa^sin top t. 

Now consider the last two terms in equation (6), the error terms. Ordi- 
narily k = c, which is small. The phase error o is also small, so that 
the first error term, which contains the product oc, can be ignored. It 
is the second error term which predominates. For w = <op, e^ = 2s/{k + q) 
cos (ot. As long as c > e, which occurs for moderately high damping, the 
value of k can be adjusted so that k + o = c, and the compensator has the 
correct output 2 ja/c cos upt at fp. For the compensator built at the Lewis 
Research Center, the value of o at w = up is about 0.1 while the typical 
drag force anemometer has a damping coefficient which does not exceed about 
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0.03. In this case c < e so k must be negative to have k + e = ;. 

This is accomplished merely by removing the inverter from the circuit of 
figure 5. The exact value for the minimum damping coefficient which can be 
accommodated without removing the inverter depends on the values of the 
extra resistor and capacitor added to the differentiator to reduce the gain 
at high frequencies and on the natural frequency of the anemometer. An 
exact formula is given in appendix B. 

The above discussion indicates that it is only necessary to adjust the 
damping control, but the more elaborate analysis of appendix B shows that 
both the amplitude and phase must be considered and that the natural fre- 
quency control must also be adjusted slightly. Nevertheless, the basic idea 
is to use a small amount of the signal from the first differentiator to can- 
cel the overall phase error which results from using the modified differen- 
tiators. 


RESULTS 

The results discussed in this section are both computed and experimental, 
and are based on the compensator circuit shown in figure 9. For the theore- 
tical curves, the settings of the natural frequency and damping controls 
were determined as explained in the example in appendix B. A spectrum analy- 
zer was used to obtain these settings for the experimental curves. The con- 
trols were adjusted by eye to give the smoothest curve on the spectrum analy- 
zer in the region of the natural frequency while the beam was excited by flow 
from an open jet with turbulence frequency components extending beyond f^. 

The first of the computed results (fig. 10) can be compared to figure 8 
to see the results of removing the inverter. The improvement in suppressing 
the large increase in gain for frequencies near f^ is apparent. 

The next two figures (figs. 11 and 12) are the computed and experimental 
results for a drag force anemometer with f^ = 1300 Hz and c = 0.024. To 
obtain figure 12, the beam was excited by the turbulence in the flow in a 
half-inch diameter free jet. 

The calculated frequency response for a high frequency anemometer is 
shown in figure 13. It has a natural frequency of 42.8 kHz and damping coef- 
ficient of 0.007. The transient response of this anemometer was also mea- 
sured by placing it in front of a shock tube, as shown in figure 14. Figure 
15 shows the output of this anemometer with no compensator in the circuit. 

The characteristic ringing is quite evident. The next figure (fig. 16) shows 
the results when the compensator is added to the circuit. The ringing is 
suppressed and the shock profile can be seen much more clearly. 

The last figure is the shock profile measured with an infinite line pres- 
sure probe, which has no characteristic natural frequency, and frequency 
response which is estimated to be flat out to the 50 to 100 kHz region. 

CONCLUDING REMARKS 

This paper has discussed the electronic frequency compensation of second 
order systems, particularly those such as the drag force anemometer which 
have small damping coefficients. 

The electronic frequency compensator has two differentiators which must 
be modified to reduce their high frequency gain. It was shown in this re- 
port that because of these modifications, lightly damped systems are not as 
adequately compensated near their natural frequencies as they might be. 
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A simple circuit modification, that of removing the inverter from the 
conventional circuit, was shown to improve the compensation near f^. . 
Numerical calculations, consisting of computer generated frequency response 
curves, and some experimental frequency response and some transient response 
curves were presented showing that adequate compensation can be achieved in 
practice. 
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APPENDIX A 
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feedback capacitor 
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ouput of first differentiator 
ouput of second differentiator 
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partial amplitude of compensator signal 
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(RfCf)-l 
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APPENDIX B 


It is shown in this appendix that it is possible to achieve acceptable 
amplitude compensation of second order system even while using differen- 
tiators modified to reduce their gain at high frequencies, so that circuit 
stability is increased and noise is reduced. What shall be done is to assume 
an input into the sensor which is to be compensated, and then examine the 
overall response of the sensor and compensator. 

It is convenient to take the input to be of the form Pq sin ut. In the 
case of the drag force anemometer, Pq is the amplitude of the variation in 
velocity head. The output x(t) of the anemometer is related to the input 
by the differential equation 

^ X X = kPq sin tot = sin ut (Bl) 


The two parameters which characterize the second order system are the 
natural frequency fp = top/2ir and the damping coefficient c. Equa- 
tion (Bl) has solution 



In terms of amplitude and phase. 
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In order to reduce the noise associated with the differentiators, they 
are modified as shown in figure 6. If ei is the input to the modified 
differentiator, and S 2 is the output, it is shown in reference 3 that 
ei and 62 are related by 

P 



(B5) 


In equation (B5), P stands for for 


r~T~> 

Rf 


and u for -s — x— 
c 


In our case, which is the input to the first 


differentiator, is the output from the drag force anemometer, x(t) = A sin 
cot + B cos cot. Then (B5) yields a differential equation for e£, which is 




CO d 

“c 


[A sin cot + B cos cot] 


= - — [A cos cot - B sin cot] (B6) 

“c 


Equation (B6) is solved by letting 62 = C sin cot + D cos cot and then 
solving for C and D, which satisfy the algebraic equations 



(B7) 


Using this solution. 
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This may be more compactly written as 



(B 9 ) 

tan e = tan [(e^ + + q>] = tan (e + <p) 

In equation (B 9 ), 


tan e, = — , tan e- = — , e = e, + 

iui rwj. i f 


The phase shift <p due to the drag force anemometer is -tan~^ B/A 
as in equation (B 4 ). 

With these definitions, 

cos 0, = , cos Qf = 

A 27 2 

wl + to/o)^ 





+ B — cos 
(0 1 
C 


cos cos C(iot - q>) - e] 



+ B — cos 0, 
"c 1 


cos sin [((i)t - <p) - (^ + e)] 


(BIO) 


The output of the first differentiator, e£, is the input of the se- 
cond, whose output is called e3. Using (BlO), e3 is easily written 
down; it is 


^3 = 




B — 2" cos ©2 cos ©^ sin [(wt — cp) — 2 {-^ ©)J 

( I ) 

c 


= - ^A^ + ^ cos^ ©2 cos^ ©^ sin [(ut - 9) - 2 ©] 


(BID 
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The output of the compensator is called e^. As shown in figure 5, if 
the inverter has been removed and k is the setting of the damping adjust- 
ment 

= ej^ + 2 k e2 + 


= Va^+ I sin (ut + 9 ) - 2 k ^ cos cos cos [(wt - 9) - 0] 

^ r 


^2 2 

cos 0jj^ cos 0^ sin (Kwt — 9) 2 0J 

^c 


= Va^ + sin (cat - 9) + T2 cos (cut - 9)] | 

where 

ta^ 2 2 

Ti = 1 - 2 k — cos cos -e^ sin -e g" cos'^ cos -e^ cos 2 e 

i ti)^ i t C I J 


(/ 2 2 

T2 = - 2 k cos Gj^ cos^ cos e cos cos g^ sin 2 g 

“c 


(B12) 




> 


(B13) 


Using the identities Va^ + B^ sin (cat - 9 ) = A sin cat + B cos tat and 

V ~~2 ? 

A + B cos (tat - 9 ) = A cos tat - B sin tat, equation (B12) can be 
written 

e^, = (ATj^ - BT2) sin tot + (BT^ + AT2) cos tat 


= Vf(^ + B^ + tI sin (tat + 9 ^) 


where 
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(B14) 


tan <pj = 


BT^ + AT2 

AT~- EJ- 


i + 


1 - 


'1 
B ^2 


-tan 9 + tan 9^. 

1 + tan 9' "tan "9^, ~ 


rhe amplitude is the product of the amplitudes due to the sensor 




and the compensator 


(yfF?) 


and the total phase shift 


9jis the sum of the two phase shifts, 9 due to the sensor and 9^, due to 


the compensator. 

Now the purpose of the compensator is to undo the effect of the probe so 
that the output of the compensator is equal to the input to the anemometer; 
that is, e£ = “• Unfortunately, perfect compensation 

with this circuit is not possible for all u. Choosing the match at w = 
however, turns out to be a good choice, since it results in excellent ampli- 
tude compensation to well beyond At u = up, then, equation (B14) 

gives 


(ATj-BT^) (BTj^AT^) 


= 0 


(B15) 


At 0) = Up, A = 0 and B = - ^ JZ c, so equations (B13) and (B15) 
give 




2 

n 2 2 

1 - 2 k — - cos cos 0^^ sin ^ cos*^ cos 0^^ cos 2 = 0 

“c 


(jj 0)^ 2 2 

- 2 k — ^ cos ©1 „ cos 0^„ cos 0„ + — T cos'^ 0. cos 0£„ sin 2 e„ = 2 c 
u_ In fn n 2 in fn n 


MB16) 


The subscript n means that the quantity is evaluated at the natural 
frequency. 

The solution of (B16) for wp/u^ and 2k is 
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(B17) 


(I) 

n 

( 1 ) 

c 




i + i: c tan 
cos cos 


sin 2 - 2 c cos 2 

cos -^1 + 2 c tan 




> 


These are the settings required for the natural frequency and damping coef- 
ficient controls to give exact compensation at u = wn. Note that if perfect 
differentiators had been used, the equations (B17) reduce to = 1 and 

2 k =s - 2 c , since then = 0. The result = 1 means that the na- 

tural frequency control is set at f,^, and the result 2 k = - 2 c means 
that the damping coefficient control is set at c. The minus sign means 
that the inverter is in the circuit. 

The second of equations (B17) provides a means of distinguishing between 
"high" and "low" damping. If the numerator is positive (c < 1/2 tan 2 en) , 
the inverter is not in the circuit and this is the low damping case. On the 
other hand, if the numerator is negative (5 > 1/2 tan 2 -On)* ^*^6 inverter is 
required and this is the high damping case. 

The value of c which establishes the boundary between the high and low 
damping cases depends on the values of the extra resistor and capacitor added 
to the differentiator circuits to reduce high frequency noise and on the na- 
tural frequency of the anemometer. This value of c is about 0.105 for 
fp =s 1300 and is about 0.18 for f^ = 42.8 kHz for the circuit shown in 
figure 9. 

The compensator schematic is shown in figure 9. It has = 505 ft, 

Cj^ = 510 pF and = 51 pF. These values were chosen in accord with 
acceptable practice (ref. 5). The value of varies from 0 to 50 kft 
to give the desired value for ug. (Fixed 250 kft resistors were added 
to allow for an fp of 1300 Hz.) 

Since the formula in (Bl7) for Wp/ug implicitly involves «p/u(- 
on the right side also, for example, through 



it is still necessary to solve for wp/ug. This can be done algebrai- 
cally but the result is a complicated cubic equation, so it is easier to 
solve by iteration. The form 



(B18) 



converges to ten places after about six iterations, starting with wo/wq = 1. 
For the example under consideration, = 2-nf^ = (2tt x 1300) 

X (505 X 510 X 10-12) = 2.104 x lO"^ and = RfCf/RfCi = CWC; = 0.1. 

For this case, o>j^l uq = 1.008. Once un/“c found, 2 k is deter 

mined from the second of equations (Bl7). For this example, 2 k = 0.1568. 
Having found wn/uc and 2 k, Ti and Jz are obtained from (B13). This is 
sufficient to calculate the overall response, as is done in figure 10. 
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Figure 1. - Drag force anemometer showing 
flow direction relative to beam: p «0: out- 
put proportional to product of velocity head 
and p. f> trf90°: output proportional to vel- 
ocity head. 













SILICON BFAM 


i C-77-2971 


Figure 2. - Drag-force anemometer. 
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Figures. ^ Blxk diagram of compensator. 
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(b) 

(a) Pure differentiator, 
(b) Modified differentiator. 


Figure 6. - Ordinary differentiator circuit and differentiator 
circuit modified for reduced noise. 
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Figure 7. - Calculated output frequency response for uncompensated 
and conventially compensated second order system (natural fre- 
quency ‘ 130Q Hz. damping coefficient - a 1) compensator controls 
set for zero output at peak of second order system. 



Figure 8. - Calculated output frequency response for uncompensated 
and conventionally compensated second order system (natural fre- 
quency ■ 1300 Hz, damping coefficient ■ a 04) compensator controls 
set for minimum output at peak of second order system. 



Figure 9. - Compensator circuit 



FREQUENCY. Hz 


Figure 10. * Calculated output frequency response for uncompensated 
and compensated (with modified design compensator) second order 
system (natural frequency - 1300 Hz. damping coefficient - 0.04) 
compensator controls set for zero output ix peak of second order 
system. 






FREQUENCY. Hz 


Figure 11. * Computed output fre^ncy response for compensated (with 
modified design compensator) and uncompensated second order system 
(natural frequency • 1300 Hz, damping coefficient - a 024) compensa- 
tor controls set for zero output at peak of second order system. 



FREQUENCY. Hz 

Figure 12. - Measured frequency response for compensated 
(with modified design compensator) and uncompensated 
second order system (natural frequency - 1300 Hz. damp- 
ing coefficient • Ol 024) compensator controls set for zero 
output at peak of second order system. 
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FREQUENCY. Hz 

Figure 13. - Computed output frequency response for compensated (with modified design compensator) and 
uncompensated second order system (natural frequency • 42. 8 kHz, damping coefficient • a 007) compen- 
sator controls set for zero output at peak of second order system. 



Figure 14. - Location of drag force anemometer in front of shock tube. 




Figure 15. -Transient response of uncompensated drag 
force anemometer to shock wave. 



Figure 16. - Transient response of compensated (with 
modified design compensator) to shock wave. 



Figure 17. - Transient response of infinite line pressure 
probe to shock wave. 
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